Ground state and excitation properties of the quantum kagome system 
ZnCu 3 (OH) 6 Cl 2 investigated by local probes. 
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We characterize the ground state and excitation spectrum of the S = 1/2, analytically pure and 
perfect kagome system ZnCu3(OH)6Cl2 using the following measurements: magnetization, muon 
spin rotation frequency shift K, transverse relaxation time , and CI nuclear spin-lattice relaxation 
T\. We found no sign of singlet formation, no long range order or spin freezing, and no sign of spin- 
Peierls transition even at temperatures as low as 60 mK. The density of states has an E 1 ^ 4 energy 
dependence with a negligible gap to excitation. 



The study of spin 1/2 quantum magnetism on the 
kagome lattice is very intriguing and lively because dif- 
ferent investigation strategies have led to fundamentally 
different predictions regarding the ground state and ex- 
citations of this system. Theories based on numerical 
or approximate diagonalization of the Hcisenberg Hamil- 
tonian favors non-magnetic (singlet) ground state [j], 
with full symmetry of the Hamiltonian Q , negligible cor- 
relation length j^j], and an upper limit on excitation gap 
of J/20 5]. These strategies provide conflicting reports 
on the possibility of a spin-Peierls state [1, Q . In con- 
trast, semiclassical treatments based on the thermal Q or 
quantum 0] "order from disorder" approach gives rise to 
a ground state with a broken symmetry of the V3 x \/3 
type, which is stable even against tunneling [lfj| . with 
gapless magnon excitations [9j, and no spin-Peierls type 
distortion [ill ]. In these circumstances, one expects ex- 
periments to provide some guidance. However, the exper- 
imental situation is equally confusing since most of the 
experimental studies of kagome-like materials suffered 
from several kinds of shortcomings, which made com- 
parison with theoretical models difficult. Some material, 
such as the jarosites [l2[ have spin S > 1/2. Others have 
5=1/2 but with a non perfect kagome structure, such 
as the volborthite Cu 3 V 2 07(OH 2 )2H 2 (CVO) A 
third class of materials such as SrCrgpGa^-gpOig (l3j 
and Ba2Sn2ZnGaio-7pCr7 p 022 [14| were hampered by 
disorder or strong third direction interaction. In light 
of these difficulties, the recent synthesis [l6[ of herbert- 
smithitc, ZnCu3(OH)eCl2, with its Cu-based quantum 
S = 1/2, analytically pure and perfect kagome lattice, 
should put the field on a new course. 

In this letter we present a comprehensive study of 
ZnCu3(OH)gCl2 using local probes. In our study, we 
address four questions which arc at the heart of the in- 
vestigation of the quantum kagome system: Do S = 1/2 
spins on kagome lattice freeze? Is the ground state mag- 
netic? What is the density of excited states, and is there 
a gap in the spin energy spectra? Finally, does the lat- 



tice distort in order to accommodate spin-Peierls state? 
We address these questions in the present work using 
nuclear magnetic resonance (NMR) and muon spin res- 
onance (/iSR) local probes. We also use magnetization 
measurements to calibrate the local probes. 

ZnCu3(OH)gCi2 was prepared by hydrothermal meth- 
ods performed at autogenous pressure. A 800 mL 
teflon liner was charged with 16.7 g Cu2(OH)2C03 (75.5 
mmol), 12.2 g of ZnCl2 (89.5 mmol), and 350 mL water, 
capped and placed into a custom-built steel hydrother- 
mal bomb under ambient room atmosphere. The tight- 
ened bomb was heated at a rate of 1 °C/min to 210 °C 
and the temperature was maintained for 48 h. The oven 
was cooled to room temperature at a rate of 0.1 °C/min. 
A light blue powder was isolated from the base of the 
liner by filtration, washed with water, and dried in air 
to afford 21.0 g product (49.0 mmol, 97.4% yield based 
on starting Cu2(OH)2C03). Magnetic and pXRD data 
were consistent with those previously reported for her- 
bertsmithite (Lsj . 

DC magnetization m measurements were performed 
on a powdered sample using a Cryogenic SQUID mag- 
netometer at temperatures T ranging from 2 to 280 K 
and fields H varying from 2 kG to 60 kG. In the inset 
of Fig. [T]we present mT/H versus T. The data collapse 
onto a single line, especially at low T, meaning that the 
susceptibility is field-independent in our range of temper- 
atures and fields. Also, no peak in the susceptibility is ob- 
served, indicating the absence of magnetic ordering. The 
only indication of interactions between spin in these mea- 
surements is the fact that mT/H decreases upon cooling 
whereas in an ideal paramagnetic system this quantity 
should be constant. Fits of the susceptibility data at 
high temperatures and lower fields than presented here 
reveal a Curie- Weiss temperature of &cw = —314 K \it\ . 
The frustration parameter Tp/\&cw\ ~ 0.22, where Tp 
~ 70 K is the temperature at which \~ l is no longer a lin- 
ear function of T, indicates strong geometric frustration. 
For comparison, in CVO T F /\@ C w \ 
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FIG. 2: A plot of the muon shift K, transverse relaxation 
time er, versus temperature. In the upper inset, the muon 
shift K versus external field H. In the lower inset, real and 

FIG. I: The muon shift K against susceptibility. In the inset, imaginary transverse field asymmetry for T = 100 K. 

normalized magnetization versus temperature. 



Muon spin rotation and relaxation (/iSR) measure- 
ments were performed at the Paul Scherrer Institute, 
Switzerland (PSI) in the GPS spectrometer with an He 
cryostat, and in the LTF spectrometer with a dilution re- 
frigerator. The measurements were carried out with the 
muon spin tilted at 45° relative to the direction of the 
beam. Positron data were collected in both the forward- 
backward (longitudinal) and the up-down and right when 
available (transverse) detectors simultaneously. Data 
were collected at temperatures ranging from 60 mK to 
200 K with a constant field of 2 kG. In GPS we used a 
resistive magnet with a calibrated field; in LTF we used 
a superconducting magnet. We took data at overlapping 
temperatures to calibrate the LTF field. The longitudi- 
nal data, which measure muon spin T± , were found to be 
useless since the T\ in our field range is much longer than 
the muon lifetime and showed no temperature variations. 

In the lower inset of Fig. [5] we show real and imagi- 
nary transverse field [TF] data taken at H = 2 kG and 
T = 100 K. The data are presented in a rotating reference 
frame (RRF) at a field of 1.9 kG. The TF asymmetry is 
best described by A TF = A cxp (-t 2 /(2T 2 * 2 )) cos(ut+4>) 
where T 2 is the transverse relaxation time, and u) is the 
frequency of the muon. The quality of the fit is repre- 
sented by the solid line. 

In Fig. [T] we depict the frequency shift, K — (uq — 
u>)/u)Q where loq is the free muon rotation frequency in 
the RRF. The difference in frequency between free and 
implanted muons is a consequence of the sample magneti- 
zation; therefore, the shift is expected to be proportional 
to the susceptibility. Indeed, as shown in the main panel 
of Fig. [TJ there is a linear relation between K and the 
susceptibility x = rn/H, with the temperature as an im- 
plicit parameter; some representative temperatures are 
shown on the upper axis. In the upper inset of Fig. [2] we 
present the field dependence of the shift at T = 10 K. 



Within the error bars, the shift in the RRF does not de- 
pend on the external field. This is in agreement with 
the susceptibility results. Therefore, Fig. [2] could serve 
as a conversion graph from muon spin frequency shift to 
susceptibility. 

In the main panel of Fig. [21 we depict K as a function 
of temperature. An additional axis is presented where 
K has been converted to \ as discussed above. We find 
that K (and hence x) increases with decreasing temper- 
atures and saturates below T ~ 200 mK at a value of 
X = 15.7(5) x 10~ 3 cm 3 /mol Cu; the error is from the 
calibration. It should also be pointed out that the en- 
ergy scale associated with spin 1/2 in a field of 2 kG is 
200 mK, and the saturation could be a consequence of 
the external field. The saturation of x is a strong ev- 
idence for the lack of impurities in our sample. More 
importantly, it indicates the lack of singlet formation or 
spin freezing. The last conclusion is also in agreement 
with neutron scattering measurements [171 ] and zero field 
HSR [H. 

The muon transverse relaxation rate 1 / T 2 * is also pre- 
sented in Fig. [2] Roughly speaking, it has the same 
temperature behavior as the shift (and as the suscepti- 
bility). T 2 relaxation is a result of defects in the sample 
causing a distribution of muons to electronic spin cou- 
pling constants or a distribution of susceptibilities. It 
has been shown that when the muon relaxation rate be- 
haves similarly to the shift 20] (or susceptibility (2lj ) 
upon cooling, it indicates quenched distribution of either 
the coupling constants or susceptibilities. In this case 
the relaxation increases simply because the average mo- 
ment size increases. Since the coupling constants and 
susceptibility are functions of distances between muon 
and electronic spin or between two electronic spins, our 
results are consistent with a lack of lattice deformation 
in ZnCu 3 (OH) 6 Cl 2 . 

Since muons could not reveal dynamic T\ information, 
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we performed 37 CI and 35 CI NMR experiments on the 
same sample. Using the two isotopes, we are able to 
determine the origin of T\. The first step in such a mea- 
surement is to find the line shape and to identify the 
isotopes and transitions. This measurement was done at 
a constant applied frequency of v avv = 28.28 MHz and 
a varying external field H . A standard spin-echo pulse 
sequence, n/2 — t — ir, was applied, and the echo sig- 
nal was integrated for each H. In Fig. [3] we show a 
field sweep for both CI isotopes obtained at T = 100 K. 
A rich spectrum is found and is emphasized using five 
x-axis and one y-axis breakers. This rich spectrum is 
a consequence of the CI having spin 3/2 for both iso- 
topes. In the case where the nuclei reside in a site with 
non cubic local environment and experience an electric 
field gradient, their spin Hamiltonian could be written as 
H = -fcK,L(l+K).H+(/ii/ )/6 [3/ 2 - I 2 + V (/ 2 - I*)] 
where vq is the quadrupole frequency, < r\ < 1 is 
the anisotropy parameter, K is the shift tensor, and 
vi = 7_ff/(27r). The powder spectrum of such nuclei has 
two satellite peaks corresponding to the 3/2 < — > 1/2 and 
—3/2 < — > —1/2 transitions, and a central line from the 
1/2 < — > —1/2 transition, which is split due to the pow- 
der average. The transition names are presented in the 
figure. The satellite peaks at T = 100 K of 35 CI are at 
6.52 and 7.07 T, and for 37 C1 at 7.91 and 8.41 T. The 
lack of singularity in the satellite spectrum indicates that 
the CI resides in a site with rj > 0, namely, with no xy 
symmetry. These findings are consistent with a 3m site 
symmetry of the CI ion in the space group R3m. 

In contrast to the two satellites, the spliting of the 
central lines at T = 100 K is clear, and appear for the 
35 C1 at 6.778 and 6.801 T and for the 37 C1 at 8.148 
and 8.161 T. Under some assumptions these values could 
be used to determine the parameter of the nuclear spin 
Hamiltonian (l3| ; assuming that the nuclear spin opera- 
tors, I x , I y and I z are colinear with the principal axes of 
the shift tensor, and that the in-plane shift is isotropic 
with K± = (K x + K y )/2, we find for both isotopes, 
K x ^ -0.0017(5), K z ~ 0.035(9) and r) = 0.4, and 
35 i/ Q = 3.75 MHz and ^ v Q = 2.55 MHz. The ratio of v Q 
is as expected from the ratio of the quadrupole moments. 
Due to the assumptions, the value of K z should only be 
considered as an order of magnitude. Nevertheless, it 
is not very different from that of the muon shift (in the 
laboratory frame) at the same temperatures. This means 
that both probes experience a similar field generated by 
the Cu spins, which for muons is usually a dipolar field. 

Temperature dependence field sweeps of the 35 CI cen- 
tral line are shown in Fig. |4j The intensities are in arbi- 
trary units for clarity. The ±1/2 < — ► Tl/2 transitions 
are easily observed at T = 300 and 100 K (indicated by 
the arrows in the figure) but are smeared out at lower 
T. In fact, the lines become so broad that the NMR 
shift cannot be followed to low temperature; hence the 
importance of the /iSR results. 

Finally, we measured the 37 CI spin-lattice relaxation 
rate T^ 1 to determine spin gap and excitation spec- 
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FIG. 3: A field sweep of 35 C1 and 37 C1. 



trum. The data were taken at a field of 8.15 T which 
corresponds to the low field peak of the central line. 
We use a saturation recovery pulse sequence. In Fig. 
[5] we depict T{~ normalized by 7 2 where 37 7 = 3.476 
MHz/T on a semi-log scale. T^ 1 increases upon cool- 
ing down to 50 K and then sharply decreases. We also 
present 35 C1 (Ti7 2 ) -1 where 35 7 = 4.172 MHz/T below 
50 K in order to determine the origin of the dynamic 
fluctuations. These measurements were done under the 
same conditions as 37 CI. When considering all temper- 
atures we find that T 35 /^ 37 = 0.75(10). From a mag- 
netic relaxation mechanism we expect this ratio to equal 
( 37 7/ 35 7) 2 = 0.69. From a quadrupole based mechanism 
we anticipate ( 37 Q/ 35 Q) 2 = 0.62 where Q is the nuclear 
quadrupole moment. Our finding is in favor of relaxation 
mediated by a magnetic mechanism as indicated by the 
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FIG. 4: 35 C1 field sweep (y = 28.28M-HV) at different tem- 
peratures. The arrows indicate the central line singularities 
observed at high-T but smeared out at low T. 
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FIG. 5: A semi-log plot of the CI inverse spin-lattice relax- 
ation, (7 2 Ti) _1 , versus temperature. Inset, a linear plot of 
the low-temperature region. The black line is a fit to Eq. [T] 
The dashed line is straight. 

overlapping (T17 2 ) -1 data points in Fig. [5] 

In the inset of Fig. [5] we zoom in on the low T data us- 
ing a linear scale. A first glance suggests that at low tem- 
perature 1 jT\ is a linear function of T as indicated by the 
dashed line. A remanent relaxation at zero temperature 
1/T™ could be due to magnetic fluctuations form other 
nuclear moments such as the protons or copper, since 
they continue to fluctuate even when the electronic mo- 
ments stop. Thus the spin lattice relaxation due to elec- 
tronic contribution only 1/Tf = 1/T X (T) - l/T[ l obeys 
1/(T£T) — const. This relation is expected in the case 
of free fermions and might be related to recent theories 

A different approach to Ti interpretation is in terms of 
magnon Raman scattering where 

1 1 f°° 

7 F {T) = 7 F ;WA 2 / p\E)-n{E)-[n{E) + l]dE (1) 



with p being the density of states, A the gap, A is a 
constant derived from the hyperfine c oup ling, and n(E) 
the Bose-Einstein occupation factor [2J|. This expres- 
sion is constructed from the population of magnons be- 
fore and after the scattering, with the associated density 
of states and the assumption that they exchanged negli- 
gible amount of energy with the nuclei since its Zeeman 
splitting is much less than a typical magnon energy. How- 
ever, in frustrated magnets the magnon might not be the 
proper description of the excitations [HI, [3] ■ Neverthe- 
less, we use Eq. Q] since it is expected for any kind of 
bosonic excitations, and since there is no other available 
theory. We assume p{E) ~ E a , with a and A as fit pa- 
rameters. The fit of Eq.[T]to the data is presented as the 
solid line in Fig. [SJ and in its inset. We find a = 0.23(1) 
and A = 0.5(2) K. Comparing to J = 209 K [ijj, this is 
a negligibly small gap. It indicates that most likely there 
is no gap in the spin energy spectra, in agreement with 
Ref. a,nd p(E) ~E X ' A . 



To conclude, susceptibility measurements down to 
60 mK suggest that there is no freezing and only a sat- 
uration of susceptibility, namely, no singlet formation. 
The data also do not support the presence of lattice de- 
formation. Finally, CI NMR T\ measurements find a 
negligibly small magnetic gap and the density of states 
p ~ E 1 / 4 . Thus, ZnCu3(OH)6Cl2 is an exotic magnet 
with no broken continuous symmetry but gapless exci- 
tations. It might be an example of algebraic spin liquid 
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